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New Cosmological Solutions in Massive Gravity 
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We find new, simple cosmological solutions with flat, open, and closed spatial geometries, contrary 
to the previous wisdom that only the open model is allowed. The metric and the Stiickelberg fields 
are given explicitly, showing nontrivial configurations of the Stiickelberg in the usual Friedmann- 
Lemaitre-Robertson- Walker coordinates. The solutions exhibit self-acceleration, while being free 
from ghost instabilities. Our solutions can accommodate inhomogeneous dust collapse represented 
by the Lemaitre-Tolman-Bondi metric as well. Thus, our results can be used not only to describe 
homogeneous and isotropic cosmology but also to study gravitational collapse in massive gravity. 



It is very intriguing to explore whether or not the 
graviton can have a mass. The first attempt to add a 
mass term to the gravity action was made by Fierz and 
Pauli [H, who considered the quadratic action for the 
graviton h^^i, in flat space with the mass term 
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The linear theory with the Fierz-Pauli mass term is 
ghost-free. However, the theory does not reproduce gen- 
eral relativity in the massless limit m — >■ 0. The extra 
three degrees of freedom in a massive spin 2 survive even 
in this limit, and therefore the prediction for light bend- 
ing is away from that of general relativity, which clearly 
contradicts solar-system tests. This is called the vDVZ 
discontinuity Q. 

As pointed out by Vainshtein Q, the discontinuity 
can in fact be cured by going beyond the linear the- 
ory. Massive gravity has a new length scale called the 
Vainshtein radius, below which the nonlinearities of the 
theory come in and the effect of the extra degrees of free- 
dom is screened safely. The Vainshtein radius becomes 
larger as m gets smaller, and thereby a smooth massless 
limit is attained. 

However, the very nonlinearities turned out to cause 
another trouble. Boulware and Deser argued that there 
appears a sixth scalar degree of freedom at nonlinear or- 
der, which has a wrong sign kinetic term, i.e., the sixth 
mode is a ghost Q. The ghost issue was emphasized in 
the effective held theory approach in Ref. [H. The pres- 
ence of the Boulware-Deser (BD) ghost has hindered us 
from constructing a consistent theory of massive gravity. 

Recently, a theoretical breakthrough in this held has 
been made. Adding higher-order self-interaction terms 
and tuning appropriately their coefficients, de Rham 
and collaborators successfully eliminated the dangerous 
scalar mode from the theory in the decoupling limit 
Then, Hassan and Rosen established a complete proof 
that the theory does not suffer from the BD ghost insta- 
bility to all orders in perturbations and away from the 
decoupling limit Q. Thus, there certainly exists a non- 



linear theory of massive gravity that is free of the BD 
ghost. 

In addition to the theoretical interests described above, 
the mystery of the accelerated expansion of the Uni- 
verse motivates massive gravity theories as a possi- 
ble alternative to dark energy. Since the attractive force 
mediated by a massive graviton is Yukawa-suppressed by 
a factor e"""", massive gravity theories with m ~ Hq 
(the present Hubble rate) could help if one were to avoid 
dark energy. Indeed, the DGP model [l^l, a concrete 
realization of massive gravity in the context of extra di- 
mensions, admits a self-accelerating solution without the 
need of dark energy 



Then, one may wonder whether or not the massive 
gravity theory developed by de Rham and collabora- 
tors [a 0] can admit flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) cosmology. Starting from the usual 
FLRW metric ansatz, it has been shown that a spatially 
flat solution is prohibited in the above massive gravity 
theory [lij. Though the same argument applies to the 
closed model as well, this interesting fact does not hold 
true for the open model, and indeed the open FLRW 
solution has been obtained in Ref. The conclusion, 
however, depends upon the form of the Stiickelberg fields 
one chooses. In other words, one can start from a non- 
standard form of the cosmological metric in the unitary 
gauge, and then move to the usual FLRW coordinates 
with a nontrivial form of the Stiickelberg fields, which 
would lead one to different conclusions. 

In this Letter, we show that fiat, closed, and open cos- 
mological solutions can indeed be realized even in mas- 
sive gravity, starting from a general Painleve-Gullstrand 
(PG) metric [3, |15| in the unitary gauge. The key 
trick used is that the FLRW metric can be recast in 



a (less familiar) PG form [l6|. Our new solutions also 
include a spherical, inhomogeneous dust collapse model 
described by the Lemaitre-Tolman-Bondi (LTB) solution 
represented in the PG-type coordinates. Thus, our so- 
lutions can accommodate not only the flat cosmological 
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model"'^ but also gravitational collapse solutions. Many 
other interesting solutions expressed in the PG coordi- 
nates in general relativity are also solutions to massive 
gravity. 

The nonlinear massive gravity theory we consider is 
described by the action 0,0] 
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where R is the Einstein- Hilbert term and U is the poten- 
tial for the graviton, 
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with two free parameters and 014, in addition to the 
graviton mass to. Each term is defined as 



U2 := 
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-6[/C2][/C]^ +8[/C3][/C] - 


f 3[/C2]2-6[/C4],(6) 


where 
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and rectangular brackets stand for traces. The tensor 
Spi, is written in terms of four Stiickelberg fields as 



^Vab, 7^6 = diag(-l, 1, 1, 1) 



(8) 



Sm denotes the action of matter which is minimally cou- 
pled to gravity. 

The equations of motion derived from the action are 
of the form 



(9) 



where X^i, represents the contribution from U due to the 
graviton mass and r^,y is the matter energy-momentum 
tensor. Here it should be noted that the effective energy 
momentum tensor X^^, from U is determined only by 
algebraic manipulation of the inverse metric matrix g^^ 
in K-i^ defined by Eq. 0. 

In this Letter, we work in a one-parameter family of 
the above theory, where 013 and 014 are given by 
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as 
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This parameter choice (made also in Ref. [l7|) enables us 
to find new solutions in massive gravity. In addition, this 
choice is useful for avoiding potential ghost instabilities 



^ Though inflation does not necessarily predict an exactly flat uni- 
verse, such a solution is useful for describing our Universe and 
investigating the cosmological perturbations. 



suggested in Ref. [l8| , as would-be dangerous fluctuation 
modes become nondynamical from the beginning. 

Our metric ansatz is taken to be the general PG 
form [iij: 

ds^ = -V^{t, r)dt^ + U^{t, r) (dr + f{t,r)dt^'^ 
-HM^^(t,r)rMrj2, (^^^^ 

where e = ±1. We write the Stiickelberg fields in the 
unitary gauge as 
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where -h is the unit radial vector, n = 
(cos (fisinO, sin ipsmO, cos 9) . 

In the absence of matter, the de Sitter solution has 
been constructed in the coordinate system of pT|) [l7| : 

ds^ = -K^dt'^ + (dr ± Hrdt^^ + a^r^dVl^, (13) 

where H = Km/\/ia, a := a/(a-f 1), and k is an integra- 
tion constant. This solution is different from the de Sitter 
solution found by Koyama, Niz, and Tasinato [l^[2^. In 
fact, the metric (|13p solves the equations of motion only 
in the special case where the parameters are given by 
Eq. (Uni). 

We are going to generalize the work of Ref. 



IJ to 



accommodate a wider class of dynamical solutions in- 
cluding cosmological ones. In light of the de Sitter solu- 
tion (|13p . we concentrate on the case satisfying 



W{t,r) 



a :- 



(14) 



The key observation here is that for any metric of the 
form with and for the Stiickelberg ilS^) . the 

tensor X^^i, reduces to the effective cosmological constant 
term: 



(15) 



Xfii, — g^iv 
a 



This happens only for the special parameter choice pi 
Consequently, any PG-type solution in general relativ- 
ity (with a cosmological constant) is also a solution to 
massive gravity. Equation (jl5|) implies that the identity 
MliTO^V^X^'' = Vf,T^"' - ^I^{s/^G^"' = is automati- 
cally satisfied. Our finding thus extends the observations 
made in Refs. 



171 . |21| to the general PG metric. 
Let us demonstrate how cosmological solutions are ob- 
tained using the above fact. The FLRW metric can be 
rewritten in a general PG form as 1^ l^j 



ds^ 



'K^dt 



1- K&^r^/a^{t) 



dr 



-rdt 



(16) 
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where K = 0, ±1, a dot denotes differentiation with re- 
spect to t, and we have taken e = — 1 for expanding flow. 
In general relativity, one is free to rescale the time co- 
ordinate to remove the constant k. However, this is not 
the case in massive gravity because such a rescaling will 
change the tensor E^i/, and then Y^^i, will depend on k. 
For this reason, we do not set k = 1 but rather leave n in 
the metric as an integration constant characterizing the 
solution. 

We include a perfect fluid whose energy momentum 
tensor is given by 



However, this coordinate transformation brings the 
Stiickclberg scalar fields to 



so that 



T 
K 



1_ _ a'H'g' 



(25) 



2 — ziy^drdg 



(26) 



= (p + p) Uf_,u^ + pgf_,u, 



(17) 



where = (—F, 0,0,0). The energy density and pres- 
sure may depend on both t and r in general, but in the 
present cosmological setting they are supposed to depend 
only on t: p = p{t),p = p{t). 

The equations of motion ([9|) now read 
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where H{t) :~ dlna/di. From the energy conservation 
equation, u^V yT^^'^ — 0, we obtain 



p + 3H{p + p) = 0. 
Rescaling the time coordinate as 



(20) 



(21) 



and using H := dlna/dr instead of H, one can appar- 
ently remove the constant k from the above equations. 
Thus, the standard cosmological equations with the ef- 
fective cosmological constant 



a 



(22) 



are reproduced. It should be emphasized that in the 
present case spatially flat, open and closed models are 
possible^ This is in sharp contrast to the findings in 



Refs. 12, l3|. It is straightforward to include a "bare" 
cosmological constant A by shifting Aoff — >■ A -I- rv? ja. 

The metric (|16p may be expressed in a more familiar 
"cosmological" form. This is done by using the new radial 
coordinate defined as 



r ^ g 



(23) 



In terms of g, the metric is indeed of the FLRW form: 

dg^ 



ds^ = -dr^ + 



1 - Kg^ 



(24) 



Thus, we see that, though the geometry described by the 
metric is spatially homogeneous and isotropic, the 
tensor T,^^ does not respect the same symmetry. After 
the coordinate transformation (|2T|) and (|23)) . J^^i, carries 
the information about the constants k and a. 

Introducing the Stiickelberg fields gives rise to a new 
invariant = g^^^ d^(lf'du(t>^ ■ It is easy to check that 
the solution we have obtained shows no singularity in 
/"^. This owes to the fact that the metric of the PG 
form has no coordinate singularity on its horizon. 

It would be an important next step to study perturba- 
tions around our cosmological background. The analysis 
of cosmological perturbations will be nontrivial because 
the reference metric does not respect the same symmetry 
as the FLRW one. However, fluctuations around the de 
Sitter background have been investigated in the decou- 
pling limit in Ref. |17| . and it has been found that the 
kinetic terms for the helicity-0 and helicity-±l modes 
vanish identically. This is the consequence of the spe- 
cial parameter choice (jlOp . In particular, as mentioned 
above, this implies that one can avoid potential ghost 
instabilities suggested in Ref. [l8| . 

Finally, we mention the description of gravitational 
collapse in massive gravity in terms of the PG-type solu- 
tions. Our cosmological metric in the PG coordinates can 
also be utilized to describe collapsing matter in spheri- 
cally symmetric spacetime (with e = +1 to adapt the con- 
tracting setup). Indeed, the metric of the PG form has 
been used to analyze the spherical contraction model of a 
star with uniformly distributed dust (22| , and with a per- 
fect fluid [l5| , while in the present case the Schwarzschild- 



de Sitter solution found in Ref. [17| is to be used to de- 
scribe the exterior solution. In addition, we explicitly 
show that the LTB metric representing general inhomo- 
geneous collapsing dust is also included in the class of 
our PG-type solutions. 

In general relativity, the metric representing spheri- 
cally symmetric spacetime can be expressed in a general 
PG form without loss of generality as 



ds" 



-N^dt^ 



52 



i + 2£; 



[Nrdt + dry + dVMO^ (27) 



where N{t,r) > is the lapse function, Nr{t,r) is the 
radial component of the shift vector, and E{t,r) > — 1. 
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The energy momentum tensor T^^ is given by Eq. p?)) 
while p and p are no longer homogeneous, and a "bare" 
cosmological constant A may also be included. The met- 
ric (P7)) can be used to analyze spherical collapse of a 
perfect fluid [24 1. 

From the equation of motion, we see 



2 V N 



M 



where wc defined an enclosed mass 

M(r, t) := Att [ {p + M|iA) r'^dr'. 



(28) 



(29) 



Here, E represents energy, which is conserved for dust, 
as shown later. 

The analytic solution for general inhomogeneous dust 
collapse is most commonly expressed in the LTB coor- 
dinates. We change the coordinates from (i, r, 6*, (j)) to 
{T,R,e,(t)) such that t = t{T) = T and r = r(T,i?) = 
f(T, R) I OL with 



-A^J^-H2i., 
r 



(30) 



where we have taken only the positive root of Eq. (pS)) for 
a collapsing fluid. Now, A^, M, and are all functions of 
T and R, and the metric is of the LTB form: 



ds^ = -AfMT^ 



f^dfl^ (31) 



In the case of dust in which one can eliminate the pressure 
gradient from the relevant equations, A^ is independent of 
R and hence is allowed to be synchronous, A^ = 1. How- 
ever, for the application to massive gravity, it would be 
better to leave A^ in the metric as an additional arbitrary 
function of T . Moreover, it turns out that E is indepen- 
dent of T, and then the remaining evolution equations 
have three types of solutions depending on the sign of 
E{R). This is the well known LTB solution that can be 
expressed in the PG form (P7)) as well. 

On the other hand, in massive gravity with the param- 
eter choice pO)) . one immediately sees that the LTB so- 
lution in the PG form ((27|) with a cosmological constant 
A + riT? ja is the solution to ^ in the unitary gauge. 
This is because, as mentioned above, given metric func- 
tions of the generalized PG form, the effective energy 
momentum tensor A"^^ of U reduces to the cosmologi- 
cal term, (I/q!)^^!/, irrespective of their coordinate de- 
pendence. The coordinate transformation from PG to 
LTB with the rescaling T T /N brings the Stiickelberg 
scalar fields to 



dT 

A^' 



-n , 



(32) 



so that 



dT2 



df dT df dR 



Note that also in this case the invariant does not 
diverge on the horizon by virtue of the PG coordinate 
system. 

In summary, we have found new cosmological solutions 
in massive gravity with flat, open, and closed spatial 
geometries. Our solutions can also describe inhomoge- 
neous gravitational collapse of dust represented by the 
LTB metric. The key was that general PG-type met- 
ric gives rise to an effective energy momentum tensor of 
a cosmological constant w? ja for the special choice of 
the parameters. This is essential to hunt for analytical 
solutions in massive gravity from the seed solutions in 
general relativity. Thus, our solutions can be used not 
only for homogeneous and isotropic cosmology with arbi- 
trary spatial curvature, but also for the spherical collapse 
model of the formation of cosmic structure such as stars 
and galaxies. 

In this Letter, the special choice of the parameters 
of the theory enabled us to have the desirable struc- 
ture cx g^y. Consequently, the conservation law 
V^X^"^ = is satisfied automatically. The presence 
of the conservation law suggests the presence of some 
accidental symmetry in the PG-type metric. It would be 
thus interesting to understand more deeply the nature 
of the metrics giving the effective cosmological term. 
As a concrete example, we will report the behavior 
of perturbations on our cosmological background in a 
separate publication [25 1 . 
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Note added While this paper was being completed, 
Ref. [2^ appeared, in which a similar situation is dis- 
cussed where an effective cosmological term X^i. ex. g^v 



arises. The authors of Ref. [26[ derived a constraint that 
is to be imposed on the metric and the Stiickelberg func- 
tions, while we have given explicit solutions for the metric 
and the Stiickelberg fields in a simple manner. 
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